Universality of equilibrium one-dimensional transport from gauge invariance 
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In this letter we address the question how interactions affect the DC conductance of a one- 
dimensional electron system not necessarily adequately described by the Luttinger model. Using 
a Laughlin type argument, we show that gauge invariance protects the universal value of the con- 
ductance of e 2 /h per channel per spin orientation if the system possesses two conserved charges 
conjugate to the chemical potentials of the external reservoirs. 



The DC conductance of a clean one-dimensional elec- 
tron system described by the Luttinger model is given 
by the universal value of e 2 /h per channel and per spin 
orientation In this letter we address the question 
whether the conductance can deviate from this universal 
value in situations where the Luttinger model may not 
be applicable. 

Such a situation may be met, for example, if one con- 
siders a system with a strong bias, so that the Luttinger 
fixed point may no longer describe the state of the sys- 
tem. In principle, one should then take into account an 
infinite number of irrelevant operators which were sup- 
pressed at the fixed point. One might expect that this 
would lead to a nonlinear dependence of the current on 
the voltage drop and to a differential conductance de- 
pending on the bias. 

Another possibility is the existence of fixed points of 
one-dimensional electron systems different from the Lut- 
tinger models. The latter arise as a result of a pertur- 
bative RG analysis for short range repulsive potentials. 
One cannot exclude the possibility that there exist other 
conformal theories describing one-dimensional electrons. 

Our analysis is based on a Laughlin type argument ||] 
and relate the quantization of the conductance in a one- 
dimensional electron system to gauge invariance. The 
only conditions imposed on the dynamics of the system 
are the existence of conserved charges conjugate to the 
chemical potentials of the reservoirs and the absence of 
impurity backscattering. Under these assumptions, we 
show that the conductance is equal to e 2 /h per channel 
and per spin orientation. 

Recently, experiments were done measuring the con- 
ductance in very pure quantum wires gj. At small 
bias, experimental data show that the value of the con- 
ductance is somewhat smaller than the universal value: 
G ~ ng{T)2e 2 /h, n = 0,1,2,..., where g(T) < 1 is a 
factor depending on the temperature T . The factor g(T) 
can be interpreted as arising from a renormalization of 
the voltage drop and must be explained in terms of the 
physics of the leads and of the junctions between the wire 
and the leads. 

To describe the conductance of a one-dimensional sys- 
tem we use an extension of the Landauer-Biittiker ap- 
proach j| to the case of interacting systems || . We con- 



sider a one-dimensional system put into contact with two 
external reservoirs with different chemical potentials, p F 
and /if/. We assume that there are two conserved charges, 
Ql and Qr, conjugate to p F and Pr, which commute 
with the Hamiltonian TL of the system and with each 
other: 



[H, Ql] = , [H, Qr}=0 ) [Ql, Qr] = 0. 



(1) 



Then the thermal state of the system connected to the 
external reservoirs is given by the density matrix 



S M = exp(-/3(H + hlQl + PrQr)) 



(2) 



and its transport properties are described by equilibrium 
statistical mechanics. 

A special property of one-dimensional electric trans- 
port is that the continuity equation for the electric cur- 
rent 



d t p + d x j = 
is solved by introducing a scalar operator 
p = ed x (f> , j = -edt<j>. 



(3) 



(4) 



Here e is the elementary electric charge. This enables 
one to derive a formula which is a one-dimensional DC 
analogue of the Kubo formula. In the equilibrium state 
characterized by the chemical potentials pl and pr the 
expectation value of the current is given by 



(j(x)) u = -e(cW(z)} M = = 
= J(TrE AI )^ 1 Tr(S M [H,^)]) = 

1 

= --^([PlQl + PrQr, (p[x)])a- 



(5) 



Formula (|5|) expresses the current in terms of the commu- 
tation relations of the conserved charges with the bosonic 
field <f> providing a solution of the continuity equation. 

Before applying formula to a general situation, we 
consider two simple examples. We start with a zero- 
coupling Luttinger model described by the Lagrangian 



C = ihip* L (d t - v F d x )ip L + ihipx(d t + v F d x )il)i 



(6) 
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The corresponding Hamiltonian looks as follows 

Ho = ihv F J dx(ip* L d x ip L - if)* R d x ip R ). (7) 
Standard bosonization formulas, 



h 
2 
h 
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ihip* L d x i(} L = ~n 2 L , 



ihip* R d x ^ R = '^n R 



(8) 



yield the following expression for the bosonized Hamilto- 
nian 



H = I dx(n 2 L +n 2 R ). 



(9) 



Coupling the system to two external reservoirs which feed 
left- and right-moving electrons into the system can be 
described by adding to the Hamiltonian an extra term 
which includes the chemical potentials n L and A*_r of the 
reservoirs: 



Hp. = Ho + ^lQl + ^rQr, 



(10) 



where Q l and Q R are conserved particle numbers of left- 
and right-moving particles: 



Q L = I dxn L , Q R — / dxn R . 



(11) 



Minimizing the expectation value of the Hamiltonian (|l 
gives 



Ml Hr 
n L = —r — , n R = -- — . 

hvp nvp 



(12) 



The total current j = ev F (n R — til) is given by the 
Landauer-Buttiker formula 



. e(/j, L - hr) _ f_ 
3 h h ' 



(13) 



where V is the voltage drop eV = hl — Mfl- 

Next, we study the effect of adding an interaction term 
described by the following irrelevant operator 



AH = sj dx(r L dM^ R dM. (14) 
The bosonized interacting Hamiltonian is given by 

Hmt =Hp+ eir 2 J dxn 2 L n 2 R . (15) 
The expectation value of Hint achieves its minimum when 

27T 2 £ 



2 Ml 
n L + t n L n R = —- — , 

hVF hVF 
27T 2 £ 2 

nR + — n R n L = -- — . 

tlVF ' hVF 



(16) 



Before computing the current in the interacting sys- 
tem we should check whether the interaction can possibly 
change the definition of the current in terms of jl and 
j R . The current is given by the general formula 



SH 



Sa(x) 



(17) 



a=0 



where a is the spatial component of the vector poten- 
tial. The latter enters the Hamiltonian through the co- 
variant derivatives iD x = iTid x + ea. For instance, the 
non-interacting Hamiltonian acquires the form 



Ho(a) =Ho + evp J dxa{ni, — n R ). 



(18) 



This implies that j — evp(nh — n R ), as we stated before. 

By applying the same procedure to the interacting 
Hamiltonian we easily arrive at 

H mt {a) = H (a)+ 

+sir 2 / dx(n\ + 2^-an L ){n 2 R - 2^-an R ). (19) 
J h h 

Therefore, the current gets an extra contribution from 
the interaction term: 

2e 

jint = ev F (n R - n L ) - e —n L n R (n R - n L ). (20) 

h 

In the ground state of the interacting Hamiltonian (p^), 
the expectation value of this expression is equal to 
e(jiL — (J' R )/h. Thus, the redefinition of the current ex- 
actly compensates the extra contribution caused by the 
interaction. 

Another example of an interacting electronic system is 
an interacting Luttinger model which is obtained from 
the non-interacting one by adding a marginal perturba- 
tion 



AH L ut = hev F I dx(iplip L )(ip R ip R ). 



(21) 



This perturbation changes some of the bosonization rules 



ihii)* L d x ij}L = 
-ifitp R d x ip R = 



h(l + e) 
h(l + e) 



(22) 



The bosonized Hamiltonian is given by 

HLut(a) = hv F J dx (j~^~( n L + n R ) + ™i n flj + 

+ev F J dxa(n L - n R ) + u L Q L + u R Q R . (23) 
The minimum at a = is achieved when 
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(1 + e)n L + en R 



PL 
hVF ' 

ft , \ , Ma. 
(1 + e)n R + en L = —- . 

hvp 



(24) 



The current is given by the same universal formula (13) 



We have seen that the interaction does not affect the 
conductance for two different types of interaction. For 
the marginal interaction, the bosonization formulas read- 
just so as to compensate the effect of interactions on the 
current. For irrelevant interactions considered above, the 
conductance is not changed because the expression for 
the current density in terms of ul and n R changes. 

Next, we turn to a proof of the general statement that 
formula (|l|) does not depend on the concrete dynamics 
of the one-dimensional electron system, as long as the 
reservoirs are coupled to the system via conserved parti- 
cle numbers Ql and Qr. 

Independently of the structure of interactions, the 
charge density has the form 



P = e{tp* L ipL +iP* r ^r)- 



(25) 



We recall that the expression for the electric current is 
sensitive to the particular form of interactions. 

It is convenient to introduce one-dimensional bosoniza- 
tion formulas for Fermi fields tpL and tpR- 



Pi 



-2-irid 



ipL = e 
, ip R = e 2 ™<^ 



(26) 



The bosonic fields <Pl and 4> R satisfy the following com- 
mutation relations 



[<f>R(x), <f>R(y)] = -^ £ ( x - y)> 



(27) 



[4>L(x),(j>R(y)] 



4tt' 



where e(x — y) = I , x > y; e(x — y) = — 1, x < y. The 
densities of left- and right-moving particles acquire the 
form 



n L = d x 4> L , n R = d x (, 



(28) 



The conserved charges Ql and Qr have the following 
commutation relations with the bosonic fields: 

[QlMx)] = ± , [Qr, My)} = ~. (29) 

One can introduce a bosonic field 4>{x) = 4>l(x) + 4>r(x) 
and rewrite the charge density as follows 



p = e(d x (f>L + d x cj> R ) = ed x cj>. 



(30) 



We note that all the commutation relations and bosoniza- 
tion rules listed above depend only on the kinematics of 
Fermi fields and are entirely independent of the dynamics 



of the model we consider. The only important assump- 
tion which we make is commutativity of the charges Ql 
and Qr with the Hamiltonian of the interacting system. 

In order to use the definition of the electric current, 
we should specify how the vector potential a enters the 
Hamiltonian. This is accomplished by a simple substitu- 
tion 



4>l{x) 
4>r{x) 



4>1(X) : 

4> a R (x) 



H{x) + Ti 

■ $r{x) - - 



dya(y); 
dya{y). 



(31) 



This rule is equivalent to replacing partial derivatives d x 
of fermionic operators by covariant derivatives. One can 
easily check it using ( p6| ) . Now we are ready to compute 
the electric current: 



j(x) = - 



Sa(x) 



= i^[H(4> L ,((> R ),<t) L (x) +<j) R (x)] 



-ed r 



(32) 



Here we used the Leibniz rule for derivatives and commu- 
tators and the following formulas obtained by combining 
© and © 



^l = -i|[^( a; ) ) ^(y) + fl ( 2/ )] ; 



Wh(x) 
Sa(y) 



-i-[<b a R (x),My) + My)}- 



(33) 



The electric current in bosonized form automatically ful- 
fills the continuity equation: 

d t p + d x j = ed t d x (j> - ed x d t (j> = 0. (34) 

can now be used to compute the 



Formulas (|29|) and 
conductance from 

= -^(\PlQl + PrQr, <t>(x)])n = 



-(Pi 



PR) 



(35) 



This confirms the universal conductance formula ( |13D in 
the general situation. 

This result can be understood by using a simple ar- 
gument similar to the one suggested by Laughlin Q 
as an explanation of the integer Quantum Hall effect. 
The most general form of the bosonized Hamiltonian 
which commutes with the conserved charges Ql and Qr 
is Ji = H.(p,r],a), where p is the charge density and 
X] = (til — nR) is the difference of the particle densities of 
left- and right-movers. Let us for simplicity assume that 
PL + PR = 0. 

The energy minimum of TL^ is achieved if 



6H 



8rj{x) 



-Apr - Mi)- 



(36) 
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The current in the system is defined by 



SH 



8a(x) 



(37) 



In order to relate the left hand side of ( |36| ) to the right 
hand side of (p7|) one applies Laughlin's argument. More 
precisely, we assume that our system is put on a ring 
which encircles a magnetic flux <f>. Let us adiabatically 
change $ by one flux quantum 



e 



(38) 



If the system was in the ground state at the beginning 
of the process it will be in a new ground state at the 
end. As the change of $ by a flux quantum can be com- 
pensated by a gauge transformation, this will be a new 
ground state of the same Hamiltonian. Formulas ( pi] ) 
imply that in this process Ql = J dxd x <f>L increases by 1 
whereas Qr — J dxd x 4>R decreases by 1. Thus, one par- 
ticle is transferred from the left-moving channel to the 
right-moving one. This implies 



~da\ p = ~dr)\ p . 



(39) 



By combining (39) with (|3 
universal result 



and (E37|) one obtains the 



SH 



Sa(x) 



2eSH 
h Srj 



-(fM L - m ). (40) 



These considerations are not applicable if the interact- 
ing Hamiltonian includes operators of the type cos(27m</>) 
responsible for backscattering Although, such oper- 
ators may be irrelevant, they break the symmetry gener- 
ated by the charges Ql and Qr, and hence formula (^) 
for the thermal state is no longer applicable. 
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